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Abstract 

One of old methods for finding exact solutions of nonlinear dif- 
ferential equations is considered. Modifications of the method are 
discussed. Application of the method is illustrated for finding exact 
solutions of the Fisher equation and nonlinear ordinary differential 
equation of the seven order. It is shown that the method is one of the 
most effective approaches for finding exact solutions of nonlinear dif- 
ferential equations. Merits and demerits of the method are discussed. 

Keywords /Nonlinear evolution equation; Traveling wave solution; Fisher 
equation; Nonlinear differential equation of the seven order 

1 Introduction 

One of the first method for finding exact solutions of nonintegrable non- 
linear partial differential equations was introduced in Ref. [1] and applied in 
Ref. [2]. However it seems to us that it was premature work. At that time 
investigators did not use the symbolic calculations and the method pQ did 
not came into notice. 
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Later we could observe appearance many other methods for finding ex- 
act solutions of nonlinear differential equations. Let us call here the fol- 
lowing methods: the tanh - expansion method [3],|3], the simplest equation 
method [S1E], the Jacobi elliptic - function method [7], the modified simplest 
equation method [814T0] , the cosh - function method [IIICG] , the Exp - func- 
tion method [331113!, the G /G - expansion method [15l[T6]. However some 
of the mentioned methods yield to the method by paper pQ . 

Currently we know four recent papers [TTH20] of application of our method 
for finding exact solutions of nonlinear ordinary differential equations. How- 
ever our opinion is that authors of these papers did not take into consideration 
all special particularities of our modified approach. 

The aim of this paper is to present our approach again and give some 
examples of application of our method. 

In fact we look for exact solution taking into account the following ex- 
pression 

N 



y(z) = J2*nQ n , (1) 



n=0 

where the function Q takes the form 

Q = T~L z' ( 2 ) 

1 + e z 

We note that the function Q is solution of equation 

Q Z = Q 2 - Q- (3) 

This equation allows us to find derivatives y z , y zz and so on. In the next 
section we present formulae for the six derivatives of solution y(z). 

We demonstrate application of our method for finding exact solutions of 
two nonlinear evolution equations. 

One of these equations is the famous Fisher equation 

u t = u xx + u(l-u) (4) 
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and another equation of the seven order in the form 



l^t ^ ^x d l^xxx ft y>xxxxx T ^xxxxxxx 0. (5) 

We use our method for obtaining exact solutions of these two equations using 
the symbolic calculations. 



2 Algorithm of our method 

The algorithm of our method has six steps. They are the following. 
The first step: reduction of nonlinear evolution equation to the ordinary 
differential equation. 

Let the partial differential equation in the polynomial form be given 

Ei(u,u t ,u x ,u tt ,u xx , . . .) = 0. (6) 

In the first step we use the reduction of the nonlinear partial differential 
equation to nonlinear ordinary differential equation taking the traveling wave 
solutions into account assuming 

u(x,t)—y(z), z = kx + ut. 

As result of this step we obtain the nonlinear ordinary differential equa- 
tion with parameters u and k in the form 

E 2 (y, uy z , ky z , u 2 y zz , k 2 y zz , . . .) = 0. (7) 

The second step: calculation of value N in formula (jT| that is the pole 
order for the general solution of equation (JTj) . 

This value is determined taking the pole order of general solution for 
Eq.(j7J). Substituting y(z) = z~ p , where p > into all monomials of Eq.flT]) 
and comparing the two or more terms with smallest powers in equation we 
find the value for N. 

We can use our method further if iV is the integer. In the case noninteger 
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N we have to use the transformation of solution y(z). For example if we 
obtain the value for N = — , where m is integer we can use the transformation 
for solution in the form y(z) = w(z) m , where w(z) is the new function. 

The third step: substitution of derivatives for function y(z) with respect 
to z and the expression for y(z into Eq.(Q. 

At this step at the beginning we substitute derivatives of the function 
y(z) with respect to z starting with higher derivative into Eq.([7|). For a 
example, if we are going to find exact solutions of the nonlinear sixth order 
ordinary differential equation we substitute at the beginning the derivative 
of the sixth order 



JV 



yzzzzzz = ^a n nQ n (Q -I) [{n 5 + 15 n 4 + 85 n 3 + 225 n 2 + 274 n+ 

n=0 

+ 120) Q 5 - (5 n 5 + 60 n 4 + 275 n 3 + 600 n 2 + 620 n + 240)Q 4 + 

+ (10n 5 + 90n 4 + 315n 3 + 540n 2 + 455n + 150)Q 3 - ^ 
-(10 n 5 + 60 n 4 + 145 n 3 + 180 n 2 + 115 n + 30) Q 2 + 
+ (5 n 5 + 15 n 4 + 20 n 3 + 15 n 2 + 6 n + 1) Q - n 5 ] 

and then we have to substitute derivatives of fifth order, fourth order and so 
on using the following formulae 



iv 



Vzzzzz = ^2 a n nQ n (Q -1) [{n A + 10n 3 + 35n 2 + 50n + 24) Q 4 + 

n=0 

-(4n 4 + 30n 3 + 80n 2 + 90n + 36) Q 3 + (6 n 4 + 30n 3 + 55n 2 + 
+45n + 14) Q 2 - (4n 4 + 10n 3 + 10n 2 + 5n + 1) Q + n 4 ] , 



N 

Vzzzz = ^a n nQ n {Q-l) [{n 3 + 6n 2 + lln + 6) Q 3 - 

n=0 

■(3n 3 + 12n 2 + 15n + 6) Q 2 + (3n 3 + 6n 2 + An + 1) Q - n 3 ] 



(9) 



(10) 
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N 

Vzzz = ^a n n Q n (Q - 1) [(n 2 + 3n + 2)Q 2 - 
-(2n 2 + 3n + 1) Q + n 2 ] , 

N 

y zz = ^2 a n nQ n (Q - 1) [{n + l)Q-n), (12) 

n=0 

AT 

t/^^a^Q^Q-l). (13) 

n=0 

At the end of this procedure we substitute the presentation of solution 
y(z) in the form 

N 

y = Y^a n Q\ (14) 

n=0 

As result of the third step we obtain the equation that has the function 
Q, coefficients a n , where (n = 0, 1, . . . , N) and parameters of co, k of Eq.([7]). 

Remark 2.1. Some authors [TTH2C)] prefer to substitute y(z) at the begin- 
ning. In this case we have to use derivatives function Q(z) as at application 
of the simplest equation method for finding exact solutions. 

Remark 2.2. If we do not know how we can find the order of pole for the 
general solution we can use our algorithm as well but in this case we have to 
to take the large value of integer N. 

Remark 2.3. If Eq.([7|) does not have some derivatives we do not substitute 
expressions for them into equation. 

Remark 2.4- If we look for exact solutions of the differential equation with 
more then six order we have to calculate additional derivatives for function 
y{z) using formulae ([1]) and 

The fourth step: finding of algebraic equations for coefficients a n and for 
parameters u, k and others of Eq. ([7]). 

In the fourth step of our method we transform the problem of finding 
exact solution of nonlinear differential equation Eq.© into the problem of 
looking for solution of the system of algebraic equations. 

Equating expressions at the different powers of Q to zero we obtain the 
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system of algebraic equations in the form 



)=0, 



(n = 0, . . . 



N). 



(15) 



TTie fifth step: solution of the system of algebraic equations. 

Solving the system of algebraic equations we obtain values of coefficients 
on, ojv-i, oo and ralations for parameters of Eq.Q. As result of solution 
we obtain exact solutions of Eq.(!7jl in the form (j!4j) . 

The sixth step: the presentation of solution y(z) of Eq.($7§ in more con- 
venient form and checking-up of solutions. 

Unfortunately there are many mistakes in finding exact solutions of non- 
linear differential equations and we have to give a good advice: we need to 
test exact solutions of nonlinear differential equations. The problem is that 
sometimes it is difficult to check-up exact solutions of nonlinear differential 
equations. We believe that in these cases we have to point out which solutions 
were checked up. 



3 Application of the method to the Fisher 
equation 



Consider the application of our method for looking exact solutions of the 
Fisher equation 



1. As result of the first step we obtain the nonlinear ordinary differential 
equation in the form 



2. In the second step we find N = 2. 

3. In the third step we substitute the second and the first derivatives of 
function y(z) into Eq. (IT7|) . In this case these derivatives can be written as 




(16) 



S k 2 y zz - u y z + y (1 - y) = 0. 



(17) 
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Vzz = Q (Q - 1) (6 a 2 Q 2 + 2 ai Q - 4 a 2 Q - oj , 



(18) 



y z = Q(Q-l)(2a 2 Q + ai ). (19) 

We have to substitute into Eq. (fIT]) the expression for y(z) as well in the 
form 

y = a + ai Q + a 2 Q 2 . (20) 
As result of the third step we have the following equation 

(6 5 k 2 a 2 - a 2 2 ) Q 4 + (2 5 k 2 a x - 10 S k 2 a 2 - 2 u a 2 - 2 ai a 2 ) Q 3 + 
+ (2 ijj a 2 — 3 S k 2 ai +A5 k 2 a 2 — 2 a a 2 — to a\ — a 2 + a 2 ) Q 2 + (21) 

+ (lo a\ + ai + 5 k 2 a\ — 2 a a{j Q + a — a 2 = 

4. Now we have to remember that Q is the function and equation has solution 
when the system is satisfied. It is possible when the system of the following 
algebraic equations can be solved 

6 S k 2 a 2 — a 2 2 = 0, 
2 5 k 2 a\ — 10 5 k 2 a 2 — 2 tu a 2 — 2 a\ a 2 = 0, 
2 to a 2 — 3 S k 2 a\ + 45 k 2 a 2 — 2 ao a 2 — tu a\ — a\ 2 + a 2 = 0, (22) 
cj ai + ai + 5 /c 2 ai — 2 a a-i = 0, 
— ^o 2 = 0. 

5. Solving the system of equations (I2"2"j) we have coefficients a 2 , a\ and a 
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in the form 

a 2 = 6 5 k 2 , ax = —65 k 2 , 

5 



1 3a; 5k 2 u 2 



a = - + — — h 



2 5 2 50 8 k 2 
and the following values u and k 



^(1,2) = ± ~7= ^(3,4) = =F 



/OS ^ ' 

Taking the formulae into account we obtain the value for u 



= -, W( 2 ) = --• 
6. Exact solutions of the Fisher equation takes the form 

1 



1 



1/(5,6) (*) = ( 1+exp( , + , o) ) . Z = k ilA Z + W(2) t 



2/(3,4) («) = 7- 7 — : r ( 2 ~ 77 7 — , \ J ' 2 = k (3A) X + UJ ( 

1 + exp (z + z ) V 1 + exp (2 + z ) / 

We check that all these solutions satisfy to the Fisher equation ( IT?]) . 
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4 Application of the method for the seven 
order nonlinear differential equation 

Let us demonstrate the application of our method for finding exact solu- 
tions of nonlinear partial differential equation of the seven order 

H C "xra ft ^xxxxx T ^xxxxxxx 0. (30) 

Taking the travelling wave solutions in Eq.f l30p we have the nonlinear 
ordinary differential equation after integration in the form 

k 

7 k 1 Vzzzzzz - ftk 5 y zzzz + a k 2 y zz + - y 3 + u y - C x = 0. (31) 

Note that if y(z) is solution of Eq.( l3"Tj) then —y(z) is also solution of this 
equation. It follows from the symmetry of Eq. fl3~Tj) . Substituting y = a z~ p 
into all terms of equation and comparing terms we obtain p = 3 and as 
consequence we have to take N = 3 in solution (1141) . 

As result of substitution of formulae (fl4"j) and ffl~3~]) - ffTUl) into Eq. fl31~|) we 
obtain the system of algebraic equations. Solving this system we obtain the 
following relations for coefficients a 3 , a 2 , a\, and for parameters of Eq. fl3ip 

a 3 = 24 ^105 7 k 3 , a 2 = -36 ^105 7 k 3 , 

(32) 

01 = 12^/105 7 ^ gQ — , a = 

83 y/j 83 y/j 

In formula (13"2"|) the coefficient 03 has signs plus and minus but we use only 



one sign. 



UJ 



a = 21 7 tf + *V-™l£, (33) 
r 83 ' 6889 7 

1343 /3 3 fc 2310 fc 3 /3 2 273 



571787 7 2 6889 7 83 



/3fc 5 -207fc 7 (34) 
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P 1 VSSOJlP (-293 + 3 iV2399) 



83 7' ^ 1660 7 



fc(5,6) = ± 



! J-830 7/3 (293 + 3z v / 2l399) 



1660 7 



eSOO^VlOS/? 3 1638A; 6 /3 2 VT05 
1_ 571787 7 3 / 2 + 6889^7 + 



120 n , 8 /ttt? — 8058A; 2 V105/3 4 

+ A; 8 a/1057- " ^r- 

83 M V [ 47458321 7 5 / 2 

So we have exact solutions of Eq.( 131j) in the form 



(35) 



(36) 



-36 V105 7A; 3 Q(z) 2 + 24 ^105 7 fc a Q(z) 3 , Q(z) 



1 + exp [z + z ) 

(37) 

where values a; and k are determined by formulae (1341) and (|35p . These 
solutions exist at relation for parameter a that is found from expression 
fl33|) . Exact solutions fl37|) at = fc(i,2) were checked-up by us. 



5 Merits and demerits of the method 

Let us discuss merits and demerits of our method. The first merit is that 
the method is simple in its application. The second merit of the method is 
that we can use the united formulae for all nonlinear differential equations 
in the polynomial form. What is more we can not calculate the value N for 
solution. However in this case it is better to take the big value of integer N. 
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The third merit of our approach is that this one allows us to obtain 
all solitary wave solutions and all one periodic solutions when we get the 
expansion of the general solution of nonlinear differential equation in the 
Laurent series [24J. 

The fourth merit of our approach is we can have solution taking the 
different functions into account. One can see that the function Q(z) can 
have different forms. For a example we obtain that using the equality 



Q(z) 



1 



1 



1 + e z 1 + e 2 



1 
2 



1 — tanh 



z 
27 J 



(3* 



We see that exact solution for nonlinear differential equation can be pre- 
sented using the hyperbolic tangents. 

Taking the autonomous form of nonlinear differential equation we can use 
the additional constant z Q in function Q(z). In this case we have to write the 
function Q(z) in the form 



Q(z) 



l _|_ e (z+z ) 



1 _ tanh [ - + ^ 
2 2 



(39) 



1 L . (Z Z\ 

- 1 - coth - + — 

2 L V2 2 



where zq is arbitrary constant and z\ is the constant as well that is 



l TC 



We see that using Q(z) we can present solution via the coth(z) functions. One 
can suggest many other trigonometric functions for description of solutions 
of nonlinear differential equations that can be obtained from function Q(z). 
Assuming 

Q(z) 



a_ 

G 



we obtain from formula (1141) expressions for the ^ - expansion method. 

It is important to note that our method gives all results that can be found 
by the Exp - function method. It is known (see,for a example Refs. [21-23J) 
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how many mistakes gave the Exp-function method. This method does not 
allow us to look for exact solutions for higher order nonlinear differential 
equations. 

Let us shortly discuss demerits of our approach. The first demerit of 
our approach is that this one does not allow us to find all solutions for the 
equation for cases when there are two and more brunches of expansion for the 
general solution in the Laurent series. However we have to note most of other 
methods do not allow to find such solutions as well. In these cases we have 
to use more difficult methods that were introduced recently in works [25H2T] . 

The second demerit of our method is that the approach does not allow us 
to look for two periodic solutions. For this aim we have to use more difficult 
methods that were developed in our recent works [251427] . 
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